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During this reporting period we have continued to investigate the 
instability discussed in Owocki and Rybicki (1984, Ap, J., 284 . J3 7 ; 
paper I). A second paper was submitted for publication describing the 
effects of line scattering and line overlap on the instability. This 
paper upholds our earlier conclusion that hot stellar winds are highly 
unstable to the growth of short wavelength perturbations, and it is 
very likely that formation and dissapation of shocks will occur and 
will be responsible for the heating of the winds to observed X-ray 
temperatures . 

Further work has been done on the problem of wave propagation in 
hot stellar winds. Numerical studies of pulse propagation have been 
undertaken in order to understand better the precise nature of the 
instability and the nature of information flows within the wind, which 
are important for understanding the critical point conditions that 
determine the properties of such winds, especially the mass loss 
rates. We have shown in the case of pure absorption line driving that 
the inward propagation of information occurs at the sound speed and 
not at the supersonic speed described by Abbott (1980; Astrophye, 
Journ. 24Q . 1183). It is not yet clear what the effective speed of 
information will be under conditions of pure line scattering. 


Publications ; 

Owocki, S.P., and Rybicki, G.B. Effect of Scattering on Instabilities 
in Line-Driven Stellar Winds. Proceedings of the 
Colloquium/Workshop on The Origin of Non-Radiative 
Heating/Momentum in Hot Stars, held at Goddard SFC, 5-7 June 
1984. NASA Conference Publication 2358 (1985). 


Papers e*. omitted: 

Owocki, S.P., and Rybicki, G.B. Instabilities in Line-Drive Stellar 
Winds: II. Effects of Scattering. Submitted to Astrophysical 

Journal, 25 March 1985. 


Papers in preparation: 

Rybicki, G.B., and Owocki, S.P. Wave Propagation in Stellar Winds 

Driven by Pure Line Absorption. To be submitted to Astrophy eical 
Journal. 


Presentations at Meetings and Conferences: 

Rybicki, G.B., and Owocki, S.P. Effect of Scattering on the 

Instability of Radiation-Driven Stellar Winds, presented at the 
165th Meeting of the AAS, Tucson, Arizona. 13—16 January 1985. 
Abstract published in BAAS, _16_, 993 (1984). 


INSTABILITIES IN LINE-DRIVEN STELLAR WINDS 


II. EFFECT OF SCATTERING 
S.P. Owocki 

Center for Astrophysics and Space Sciences, 
University of California at San Diego 

G.B. Rybicki 

Harvard-Smithsonian Center for Astrophysics 


Received: 


Submitted to The Astrophysical Journal 



ABSTRACT 


We extend our earlier analysis (Owocki and Rybicki 1984; paper I) of the linear 
instability of line-driven stellar winds to take proper account of the dynamical effect 
of scattered radiation. The principal differences from the previous study of scat- 
tering by Lucy (1984) are: 1) We use the physically more reasonable assumptions 
of complete redistribution (vs. coherent scattering) and spherically symmetric (vs. 
plane parallel) flow in the wind. 2) In addition to examining Lucy’s "drag effect” 
of the mean scattered radiation field on velocity perturbations, we also investigate 
the role of perturbations in the scattered field itself. 3) We quantitatively examine 
the effects of line overlap. 

Our principal findings are: 1) The drag effect of the mean scattered radiation 
does indeed greatly reduce the contribution of scattering lines to the instability at 
the very base of the wind, but the instability growth rate associated with such lines 
rapidly increases as the flow moves outward from the base, reaching more than 50% 
of the growth rate for pure absorption lines within a stellar radius of the surface, and 
eventually reaching 80% of that rate at large radii. 2) Perturbations in the scattered 
radiation field may be important for the propagation of wind disturbances, but they 
have little effect on the wind instability. 3) The contribution of strongly shadowed 
lines to the wind instability is often reduced compared to that of unshadowed lines, 
but their overall effect is not (as suggested by Lucy 1984) one of damping in the 
outer parts of the wind. 

The primary conclusion we derive from these results is thus that, even when all 
scattering effects are taken into account, the bulk of the flow in a line-driven stellar 
wind is still highly unstable. 


2 



I. INTRODUCTION 


In the first paper of this series (Owocki and Rybicki 1984; paper I), we showed 
that a flow driven by line absorption of a point source of continuum radiation is 
unstable to a velocity perturbation with a wavelength A approximately equal to or 
less than a Sobolev length L (over which the mean flow velocity increases by a ther- 
mal speed), but becomes asymptotically stable to perturbations with A » L. We 
thereby reconciled conclusions of the many earlier studies that had also found such 
flows unstable (Milne 1926; Lucy and Solomon 1970; Carlberg 1980; MacGregor, 
Hartmann, and Raymond 1980, hereafter MHR) with the work of Abbott (1980), 
who found them to be marginally stable. From this analysis, we then argued that 
the winds of hot, luminous, early-type (i.e., spectral classes O and B) stars, which 
are thought to be driven by line-absorption of the star’s continuum radiative mo- 
mentum flux (Lucy and Solomon 1970; Castor, Abbott, and Klein 1975, hereafter 
CAK; Cassinelli 1979; Abbott 1980, 1982; Cassinelli and MacGregor 1985; Abbott 
and Lucy 1985), are unstable to short wavelength perturbations. 

Sn<~h winds are actually driven mostly by strong scattering (vs. pure absorption) 
lines (Abbott 1980; Abbott and Lucy 1984), but, in the supersonic portions of the 
wind, the scattered radiation can be shown (Sobolev 1957, 1960; Castor 1970, 1974; 
Lucy 1971) to have negligible dynamical effect on the mean flow in comparison 
with the direct “core” radiation from the stellar photosphere. On this basis, we 
argued in paper I that the scattered radiation field should have negligible dynamical 
effect on the perturbed flow as well, and hence that an analysis based solely on 
the extinction (i.e., neglecting scattering) of the core radiation would suffice to 
study the stability of the wind. Subsequently, Lucy (1984) criticized this neglect of 
scattering and showed that the mean scattered radiation field exerts a liag force on 
velocity perturbations which counteracts the amplifying effect of the force associated 
with the direct extinction of the core radiation. Under the specialized assumption 
that the scattering is coherent, Lucy found that this drag effect exactly cancels the 
contribution of isolated, pure scattering lines to the instability near the wind base, 
where the flow is nearly plane-parallel. He further conjectured that “in the outer 
parts of winds, line drag results in a damping contribution from resonance lines 
that are strongly shadowed by neighboring resonance transitions”, and he thereby 
concluded that “line-driven winds would seem therefore not to be as highly unstable 
as hitherto claimed”. 

In order to examine this question of the effect of scattering lines on the stability 
of radiatively driven stellar winds, we have extended the analysis of paper I to take 
proper account of the dynamical effects of scattered radiation. Our analysis (§11 
and § III) differs from Lucy’s in that it is based on the physically more reasonable 
assumptions that the scattering occurs with complete redistribution (vs. complete 
coherence) and that the mean wind flow is spherically symmetric (vs. plane paral- 
lel); further, in addition to Lucy’s drag effect of the mean scattered field on velocity 
perturoations, we include the effects of perturbations in the scattered field itself. We 
also examine the effects of photospheric limb darkening (§IVa), photon destruction 
(§IVb), and line overlap (i.e., line shadowing; §IVc). The results (summarized in 


3 



§V) of this more complete scattering analysis lead to markedly different conclusions 
about the importance of scattering effects in reducing the instability of radiatively 
driven stellar winds. In particular, our primary conclusion is that, even when all 
line scattering effects (including line drag and line overlap) are taken into account, 
the bulk of the wind flow is still very unstable, a conclusion that is in agreement 
with that of paper l and in contradiction with that of Lucy (1984). 


II. PERTURBATION ANALYSIS FOR A SCATTERING LINE 

a) Assumptions 

In examining line-driven stellar winds for stability against perturbations, it is 
important to calculate correctly the perturbed value of the line force. In this section, 
we calculate this perturbed line-force under the following assumptions: 

1) The iine is a pure scattering line so that the photon destruction probability per 
scattering 6 = 0. 

2) The line is isolated so that its line center frequency is displaced from its near- 
est (blueward) neighbor by a velocity separation Av > 2uoo, where is the 
asymptotic wind speed. 

3) The line is optically thick, but unsaturated in the mean flow, so that scattering 
occurs in the Doppler core of the Voigt profile. We therefore ignore the line 
wings. 

4) The line scattering is isotropic. 

5) The line scattering occurs with complete redistribution in frequency. 

6) The radiation from the underlying stellar photosphere is not limb darkened. 

7) The mean wind flow is radially directed and spherically symmetric. 

8) The wind is highly supersonic so that the mean flow speed vo » Vth (where 
vth is the ion thermal speed), implying that the Sobolev approximation, L = 
w t/i/ w o ^ vol v 'o = H (« R*, the stellar radius), is satisfied for the mean flow. 

9) Perturbations in the line opacity k can be ignored. 

10) The perturbations are of the radial, sinusoidal form exp(t(fcr — w£)) with fre- 
quency u) and wavenumber k. 

11) The perturbation wavelength A i?* so that a WKB analysis is justified. 

12) The perturbations are small amplitude so that a linear analysi. is justified. 

The effect of relaxing assumption (3) has already been examined in paper I, 
where we showed that unsaturated optically thick lines make the dominant con- 
tribution to the flow instability. Relaxation of assumptions (l), (2), and (6) is 
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discussed in §IV below. Subsequent papers in this series will examine relaxation of 
various remaining assumptions, particularly (12). 


b) The Transfer Equation in the Co-Moving Frame 

We begin by considering the transfer equation in a co-moving frame represen- 
tation for spherically-symmetric outflow (assumption [7)), in which we include the 
effect! of sphericity in the frequency derivative, but not in the spatial gradient (See 
Rybicki 1970 for a discussion of this approximation.), 


81 _8I ./ \ f r cl 

m - t - - = -x4>{x) [/ - 5 ] . 


dr 


dx 


Here 


g = — 

Vth 


/xV + (l - n 2 ) ^ j 


(1) 


( 2 ) 


r is the radius, cos“ V is the angle betv’een the ray and the outward radial direction, 
x is the frequency displacement from line center measured in the comoving frame 
in units of the Doppler width u is the flow speed, v' = du/dr, <f>(x) is the line 

profile function, I is the specific intensity, and the line strength x is related to the 
line opacity k through the mass density p by x — Pl- 
under assumptions (1), (4), and (5) of pure, isotropic scattering with complete 
redistribution, the line source function S is given by the mean line intensity, 

s = 7 = (7), (3) 

where the bar denotes an average over the line profile 


_ f 00 

/(r,/i)=/ dx<f>(x)I{r>n,x), 

J — oo 

and the angle brackets denote an average over angle, 

( 7 ) s 2 J_ i d f i7{r,n). 

Finally, the line flux, 

H = {«/) . 

determines the line force per unit mass g through 


9 = 


V 


V 


H. 


(4) 


(5) 


(6) 


(7) 


i) The mean line force and scattered radiation field 

From the above we see that, to obtain the response of the line force to pertur- 
bations, we must obtain the perturbed line flux. Toward this goal, let us first divide 
the flow into a zero order mean state plus a first order sinusoidal perturbation, 


w {r ,0 = v 0 (r) + 


( 8 ) 
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with corresponding mean and perturbed values of other variables similarly denoted 
by w 0” and “1” subscripts, respectively. The zero order transfer equation is then 
given by 

Qo- 0 ^ = Xo<f>( x ) [Jo ” S 0 j , (9) 

where, in accordance with the Sobolev approximation (assumption (8]), all spatial 
derivatives of mean flow variables are ignored in favor of frequency derivative terms 
(Sobolev 1957, 1960), that is, ( dIofdr) x — 0. The Sobolev optical thickness in 
direction y is defined by 

1 _ Qo _ 1 ~ M 2 + £ 

Ty X0 T 0 Tj ’ 

where 1 r 0 = XoTVth/vo and T\ = XoVth/vb are the Sobolev optical depths in the 
transverse and radial directions. Alternatively, defining (Castor 1974), 



we may write 


TM _ d lnv 0 (r) 
° {r) - dlnr h 

in) 

to 

Tfl 1 + oy 2 ‘ 

(12) 


For an isolated line (assumption [2]; The effect of line overlap is treated in 

§IVc.), the intensity on the blue wing of the line is fixed by the intensity from the 

underlying stellar core, 

7(r,/i,x -*■ oo) = I c D(y,r), (13) 

where I c is a constant specifying the magnitude of the core intensity and D(y, r) is 
a function specifying its angle dependence at radius r. Under assumption (6) of no 
limb darkening (cf. §IVa), D(y) is unity for y > y* and zero for y < y., where 

M r ) = l/ 1 “ (^r) . ( 14 ) 

and where R* is the radius of the photosphere. 

Using the integrating factor exp(— r M $(z)), where 

$(x) = J <f>(x')dx\ (15) 

we may integrate (9) with the boundary condition (13) to obtain 

Jo = J c D(y)e~ r ^ x '> + S 0 (l- . (16) 

1 An exception to our notational convention, the subscripts on To and rj here 
refer to y = 0 and y = 1, and not to the mean and perturbed states. 
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The frequency averaged intens ity in the mean flow is 


Jo = I c D(fi)p(fj.) + S 0 (1 - p(/z)) , (17) 

where 

p[fi) = f dx<f>(x)e~ T ^^ = (18) 

J-oo Tp 

is the angle-dependent escape probability. The line flux in the mean flow is thus 

H 0 = (p7 0 ) = h {pD{p)p(p)) , (19) 

where the terms proportional to Sq have vanished because they are odd functions of 
p, with the result that the mean line force go ~ Ho (see equation (7)) is independent 
of the scattered radiation field (Sobolev 1957, 1960; Lucy 1971; Castor 1974). Under 
assumption (3) of a strong line with » 1, equation (19) becomes, 


H o 



where 

J_ _ 1 + g(l + lij) /2 _ l f J_\ 
T* ” To ~ 2 \Ti + 7>. J 

Using equation (7), the mean line force is thus 


9o = w I e 


/^*\ 2 kAuq 
\ r / cr* 


( 20 ) 


(21) 


( 22 ) 


We now turn to the problem of solving for the mean scattered radiation field. 
From equation (3), the mean source function is 


& 

II 

o 

II 

(23) 

where 


0 = <p(m)> 

(24) 

is the angle- aver aged escape probability and 


0 c = (D(p)p{p)) 

(25) 


is the core penetration probability. For the case of interest here of a strong line 
with » 1, we find, 


So Pc _ ( 1 ~ P« 2 \ / 1 + 

I c ~ 0 ~{ 2 ) \l + a/3) ’ 


(26) 


7 


where we have introduced the notation, 

= (27) 

n 1 - /u. n j”j 

Finally, we write the solution to the zero order transfer equation in the form, 


Io - S 0 = I e e~ r ^ {x) 


D -J 


{28) 


ii) The perturbed line force and scattered radi'niion field 

Next we consider the problem of solving for the perturbed scattered radiation 
field. The first order transfer ec nation for the perturbed intensity Ii is 


ikp.Ii - Qo-^f ~ Qi^~ 7 = ”Xo <f>(z) \h - Si] , 


(29) 


where we have used assumption (9) and the WKB approximation (cf. paper I) in 
ignoring a term proportional to xi« In addition, 


< 3 . 


Vth 


ikfx + 


2 , 1 “ M 2 


» ikfi 


2 ^1 
vth' 


(30) 


the latter approximation also being valid in the WKB limit fcr > 1. Using equations 
(9), (10), and (28), we write equation (29) as, 


+ (V + »>/r„/i) = r„rf(*) 5 , + (d - j ) ^ 


where we have defined 


Kfi ~ Q oifiY 


(31) 


(32) 


Then with the integrating factor exp(-tpKpX + r>^( x )) and the hlue-wing bound- 
ary condition /, (r, fi, x — * oo) = 0, equation (29) can be integrated to yield 

J,(r,At,z) = Sir^*'**-^*) f^° dx , 4>(x , )e- i ^ lt,+r ^^ + 

+I c tQ± e WM*-**(*) dx'^e-^^’ . (33) 

Averaging over frequency then yields 

7i (r, ft) = SiA + I e 


D ~J 


9ln 

Q o B ' 


(34) 
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where 


A(n) ~r h f°° [°° dx , <j>(x , )c-' ttl< '‘ x,+T »* (x, '> 

7-00 J z 


and 


B{fj.) = Tfi r dx4>{x)t iilK ^ x - r ^^ f°° dx 1 tf>(x , )e~ iftK, ‘ x ' . 

7 — 00 7 1 

The numerical evaluation of A{n) and £(/z) ia discussed in the Appendix. 
Averaging equation (34) over angle, we obtain 

7, -s. (*>+/.((*- £)af). 


(35) 

(36) 


(37) 


Then, since by equation (3), S\ — Jj, equation (37) may be solved for the perturbed 
source function, 


Si = 


h 

1-(A) 




(38) 


which may be used to eliminate S\ from equation (34). Finally, taking the first 
angle moment of equation (34), we obtain the perturbed flux, 


H x 



QiA) \ QiB\ 

i -<A)J Qo /' 


(39) 


Note from equation (35) that, since r M and Kp are even functions of /z, the real 
(imaginary) part of A is an even (odd) function of /z^JThus, {A) is purely real, while 
(fxA) is purely imaginary, and the dependence of H\ on the perturbed scattered 
field can be written in terms of the real ratio, 


(fiSSA) _ —i(fiA) 
’’ = 1 - (SU> _ !-(/!)■ 


(40) 


Since g ~ H, the perturbed line force g j can be written in terms of H \ and 
mean flow quantities, _ 

9± _ Hj_ 

9o Ho 

From equations (20), (30), (32) and (39), we then find 



£i 

Vl 



(m + in) . 


(42) 


The evaluation of the required angular averages is discussed in the Appendix. An 
approximation to g\/vx can be obtained by using equation (A3), yielding 


s-vsru”- 


2 r\ 2 /( 0 C \ M 2 [SKpIfr.) - r?] + t> 3 [l + t}Kp/{ 2z.)|' 


& 


1 + [fJ.Kpl{2x.)\‘ 


(43) 
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where u>. = 2x.go/v t h and A. = x* _l = r./( 2x.xo) are generalization* of the 
growth rate and bridging length of paper I, now derived for the case that t ho stellar 
core has a finite angular size ( /r. ^ 1). As shown in paper I, the real parts of the 
perturbed force to velocity ratio, 5i(0i/vi), essentially 2 gives the contribution of the 
line to the growth rate of flow perturbations; therefore, wo shall henceforth refer to 
the ratio simply as the “line growth rate." 


III. RESULTS 


To evaluate the perturbed force at a specific radius r in the wind, we must 
assume a specific velocity law in order to determine terms which depend on a == 
[d In vo fd in r) - 1. We may then use straightforward numerical angle quadrature to 
calculate from equation (43) the perturbed line force as a function of wavenumber 
k. For all detailed calculations presented here we assume a velocity law of the form 


vo(r) = Von 




M 


for which 


For a velocity power index p = the outward acceleration of the wind is every- 
where proportional to the local inward acceleration of gravity, in agreement with the 
theoretical result of CAK; however, observations appear to suggest (Abbott 1982) 
a less steep velocity law with p « 1-2. We Bhall show that any of these choices for 
p yield qualitatively similar results for the wind instability. 


+ l)/2. 


( 45 ) 


a) Optically Thin Perturbations 


In the very short wavelength, optically thin perturbation case, k » Xi we expect 
that the perturbed radiation field will be very small due to spatial averaging over 
the perturbed source function S\. This is confirmed in the Appendix, where it is 
Bhown that t) -* 0 in this limit. Using this fact and equation (A17), we see from 
equation (42) that the perturbed force becomes predominantly real and independent 
of scattered field perturbations. Using equations (20) and (41), we can express the 
perturbed force ratio in the form, 


£i 

Vl 


9? 


(5r)-(M)£)/(^) 

(l + < 7 / 3) (J ?3 + oEs) — (1 + oEz) (1/3 t «7/5) 


= U>' 


1 + /i. 


(1 + 0 / 3 ) (1 + a (1 + / i 2 ) / 2 ) 


{ 46 ) 


with £3 and £5 defined in equation (27). The first term in the numerator results 
from the direct extinction of “core” radiation from the stellar photosphere, while 
the second term results from the drag effect of the mean scattered radiation field. 

2 I.e., exactly for the case of zero sound speed and within a factor of 2 otherwise. 
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Figure 1 shows the radial variation of /vi ) as calculated from equation (46) 
for the velocity laws given by equation (44) with p — j, 1 and 2. Comparison with 
the analogous curves for a pure absorption line shows how the scattering-line growth 
rate, as given by 9?{ffi/vi) (see paper I), is reduced by the drag effect of the mean 
scattered radiation field So- Near the base of the wind, the growth rate is essentially 
zero; but as the flow moves outward from the base, it rapidly increases, reaching 
50% of the rate for pure absorption lines within a stellar radius from the base, and 
eventually approaching 80% of that rate far from the star. This still implies that 
perturbations in a wind driven by scattering lines will grow by many e-folds within 
a characteristic wind outflow time. We conclude that the bulk of such a wind is 
still very unstable to short wavelength, optically thin perturbations, even when the 
drag effect of the mean scattered radiation field is taken into account. 

b) Dependence on Perturbation Wavelength 

We now wish to establish the behavior of the perturbed force law g i/vj as 
a function of wavenumber k , thus establishing the “bridging law" for the case of 
scattering lines. 

As the first step, we evaluate by numerical quadrature the quantity q, defined 
in equation (40), using equations (A9) and (A10) of the Appendix. Results as 
a function of normalized wavenumber fc/x* are given in figure 2 for a line with 
radial optical depth T\ = 100 at various radii r in a wind with a velocity law 
given by equation (44) (with p = ^). Note that r? — » 0 for both k x» an ^ 
k 'JP x.ti ss x* in accord with our physical expectation that the effect of the 
perturbed scattered radiation field vanishes in both the Sobolev and optically-thin 
perturbation limits. (See also the discussion in the Appendix.) The perturbed field 
can therefore, at most, affect the instability of perturbations with a wavelength A 
smaller than the Sobolev length L (sy l/x*)> but larger than a line-center mean- 
free-path 1/x (« L/rj). Note also tnat 77 < 0 for all k. As we shall see, this means 
that the perturbed scattered field typically causes a slight increase in the strength 
of the instability, in contrast with the decrease caused by the mean scattered field. 

Using numerical quadrature to evaluate the angle averages in equation (43), 
we may now calculate the variation of the perturbed force with wavenumber k. In 
figures 3a-3d we plot for a line with radial optical depth rj = 100 the results of such 
a calculation at selected radii, r/R , = 1.01, 1.5, 2.5, and 10 in a wind with p = 5. 
The perturbed force to velocity ratio g i/tq is scaled here in units of u. while the 
wavenumber k is scaled in units of x*- For comparison, we also plot this perturbed 
force ratio in the absence of any scattered field (i.e., setting both So = — 0), as 

well as in the absence of any perturbations in an existing mean scattered field (i.e., 
setting only S 1 = 0). In each case, the solid curves denote D?(ffi/vi), and the dashed 
curves denote 5(gi/ui). Recall again (see paper I) that, while 3?(ffi/vi) essentially 
det ermines the growth rate of perturbations, 9'(pi/ui) determines the speed of their 
phase propagation. 

Concentrating first on the dashed curves in figure 3, note that the mean scat- 
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tered field Sq has no effect at all on &(<7i/t>i). This can be understood directly from 
equation (43), for which we note that the the imaginary part is odu in p, so that 
only the term representing the effect of direct radiation from the -»tar contributes 
to Q'tffi/i'i). On the other hand, perturbations in the scattered field do have a 
dramatic effect on Q(g i/t>i), reducing it to negative values for moderate and short 
wavelength perturbations (A A. « L). (Because of the logarithmic scale of the 
ordinate, we have actually plotted |&(ffi/ui)j in figure 3). As we shall see in the 
next section, this change in sign in ^(ffi/ui) results in a like change in sign in the 
phase propagation of radiative-acoustic waves of moderate to short wavelength. 

Concentrating next on the solid curves, note from figure 3a that including the 
mean scattered field S 0 substantially reduces the magnitude of 9?(^i/t/j) for short 
wavelength (A <£ A* ss L) perturbations near the base of the wind (r = 1.0172.) (cf. 
curve with So = 0 vs. those with 5b ^ 0). As discussed above, this reduction results 
from the “photon drag” exerted on the velocity perturbations by the mean scattered 
field. Comparison with figures 3b-3d again shows that this drag effect quickly 
becomes less important as the gas moves away from the wind base (r > 1.572.). In 
addition, note that including the perturbed scattered radiation field S \ or ’y slightly 
alters 92(^1 /vi), generally increasing it somewhat for long wavelength (A £ A. « L) 
perturbations near the star, but otherwise hardly changing it at all. Thus, insofar 
as the instability of the wind is concerned, we see that the overall effect of including 
scattering is to reduce the growth rate substantially very close to the star, but only 
slightly far from the star. 


c) Linear Stability- Dispersion Analysis 

We now apply this calculation of the perturbed force associated with a scattering 
line to a linear stability-dispersion analysis in which we determine the growth and 
propagation characteristics of small-amplitude perturbations in a flow driven by 
line scattering. For a small-amplitude velocity perturbation v\ with associated 
perturbed force g i, the stability-dispersion equation relating the frequency u to the 
wavenumber k can be written (see paper I) 


ui 2 — t — — qa 2 fc 2 = 0, (47) 

vi 

where 

2 w 

1 ' 3 w + iW 

describes the dependence of the perturbed pressure gradient force on the frequency 
a/ and the temperature relaxation rate W . For the sake of brevity, we will only give 
details for the isothermal ( W w; qr = l) and adiabatic ( W co; 7 = 5/3) limits. 
(The results for the intermediate damping case, u> « W, are quite similar to those 
of paper I.) Since the factor qf in both these extremes is independent of frequency, 
we may use the quadratic formula to solve equation (47), 



w± 


igi/vi ± \f^a^k 2 - (gi/t>i)‘ 


(49) 
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where g\/v\ is calculated from equation (43) using the velocity law (44). 

In figures 4a and 4b, we plot the real (solid curves) and imaginary (dashed 
curves) parts of the two frequency modes u>+ and w_ (scaled by u >. ) vs. the 
wavenumber k (scaled by x*) for various values of the sound speed ^/ 7 a (scaled by 
w./x.) at the typical stellar wind radius r = 2.5 R*. The results are quite similar to 
those plotted in figure 3 of paper I. The major difference is found in S?(w±), which 
changes sign for a wavenumber k » 2x»- (Because of the logarithmic scaling of the 
ordinate, we have actually plotted |SR(w±)|). This indicates that short wavelength 
perturbations in a flow driven by scattering lines have their phase propagation in 
the opposite direction from those in the absorption line driven case studied in paper 
I. This effect has important consequences for how information about a perturbation 
propagates through the wind, a topic which we shall examine in more detail in a 
future paper of this series. 

In contrast to this qualitative change in !K(w±), note that the variation of the 
growth rate 0(w±) with wavenumber k is only slightly reduced relative to the pure 
absorption case. The major effect of scattering on the instability of the wind is thus 
to reduce the growth rate of short wavelength perturbations as a result of the drag 
effect of the mean scattered radiation field. 


IV. DISCUSSION OF OTHER EFFECTS 


a) Effect of Limb Darkening 

We may easily generalize the above analysis of the perturbed line-scattering 
force to include the effects of limb darkening of the incoming stellar photospheric 
“core” radiation field. In this case, the angular dependence of the core intensity (cf. 
equation (13)) for fi > fx t is given by 


D{fij - 


(50) 


/(*.,!)' 

Here I is the emergent continuum intensity from the stellar photosphere as a func- 
tion of v ~ cos(0.), where 0 , is the angle between the line of sight and the outward 
normal at the stellar surface. From elementary trigonometry, we may write u as a 
function of /i as 


t/(^,r) = 


fl 2 - fj.. 


4 - 


(51) 


1 - H* 2 

With the minor modification of using equation (50) for D(fx), the analysis for the 
perturbed line force then proceeds as in §11. 

In general, the limb darkening law depends on the continuum frequency relevant 
to the particular line unaer discussion. For purposes of illustration, let us choose a 
limb darkening law to be the same as for the frequency-integrated Eddington limb 
darkening law, so that 

D(fi) = | (i -bt^r)) , (52) 
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(cf. Mihalas (1978), equation (3-37)). For this case, we find that the effect of 
limb darkening on the $(31/^1), and hence on the phase propagation of radiative- 
acoustic waves, is negligible. We further find that the effect on the JR(ffi/vi), and 
hence on the instability growth rate, is also small, except that now W(<?i/vi) « 
0.026 tu* > 0 at r = J2«, so that a wind driven by scattering lines is unstable, even 
at its base If, on the other hand, the photospheric field were limb brightened (which 
is perhaps unlikely), then flow perturbations at the base of such a wind would be 
mildly damped. 

It is interesting to note in this regard that, in the coherent scattering analysis 
of Lucy (1984), the cancellation of the instability is specific only to the assumption 
that the flow is plane parallel, and is wholly independent of any assumption about 
the angular dependence of the incorrug radiation in the blue wing of the line. 
Here we see that, under the physically more realistic assumption of scattering with 
complete redistribution, the relative importance of the drag effect depends critically 
on the angular dependence of the stellar photospheric radiation field, with exact 
cancellation of the instability occurring at the base of the wind when there is no 
limb darkening. Because of this sensitivity of the cancellation to assumptions about 
the line frequency redistribution and the angle dependence of the radiation field, 
a reliable answer t.o the question of whether the wind base is, in fact, unstable 
may require a more sophisticated treatment involving full angle- and frequency- 
dependent redistribution functions for the line scattering. This topic is currently 
under investigation. 


b) Effect of Photon Destruction 

We may also readily generalize the analysis of §11 to include the effects of photon 
destruction. (For simplicity, we still ignore the effect of any thermal creation of 
photons.) If the probability of destruction per scattering of the line photon is e, 
then the line source function relation, equation (3), becomes, 


S = (1 - e) J. 


(53) 


We may therefore still use equation (43) for the perturbed force if we just modify 
the definition of terms that arise from the mean and perturbed scattered field. We 
accomplish this by changing the definition of the core penetration probability (cf. 
equation (25)) 

Pc - (1 - e) (D(fi)P(p )) , (54) 

and the perturbed scattered field ratio (cf. equation (40)), 


»? — (! — 


X (m&A) 

} (1 - I RA> ' 


(55) 


Note then that for e = l, we obtain p c — rj = 0 and thereby recover the pure 
absorption results (cf. curves labeled So — Si = 0 in figure 3). In general, the 
results for the case of 0 < £ < 1 are intermediate between those for the c = 0 and 
e = 1 cases as plotted in figure 3. 
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c ) Effect of Line Overlap 

Many of the lines which drive stellar winds are not effectively isolated from 
neighboring lines (Abbott 1980, 1982; Olsen 1982; Friend and Castor 1983; Abbott 
and Lucy 1985), as we assumed in §11 (assumption 2); instead they often “over- 
lap” within the wind because their rest wavelengths A 0 *re separated by less than 
the maximum Doppler shift, 2Aot>oo/ c * between two points in the wi"d. In justify- 
ing his contention that line-driven stellar winds are not as unstable as previously 
thought, Lucy (1984) has asserted that in the outer parts of winds, line drag re- 
sults in a damping contribution from resonance lines that are strongly shadowed 
by neighboring resonance transitions. We can examine this assertion quantitatively 
by generalizing the perturbation analysis of §11 to the case of a pair of optically 
thick overlapping lines with a wavelength separation in velocity units of Aw, where 
0 < Aw < 2voo. (Further generalization to sets of three or more overlapping lines 
is straightforward, but tedious.) 

We first note that the above analysis for an isolated single line can be applied 
without modification to the bluer line of the pair. This is because radiation scattered 
in the wind always appears red-shifted when viewed from another location in the 
wind, and so scattering by the red line can have no effect on the blue line. On the 
other hand, for the red line we must generalize the single-line analysis to include 
the new possibility that radiation scattered by the blue line is subsequently red- 
shifted into resonance with the red line. In this case, the boundary condition for 
the intensity on the blue wing of the red line is now (cf. equation (13]) given by the 
intensity scattered out of the red wing of the blue line. Letting R = R(r, /z,Av) 
denote the blue line resonance point corresponding to the red line resonance point 
r, the boundary condition can be written 

/ r (r, — ► oo) = h{R, p\x — oo). (56) 

Defining 1 = r — R, we then have p = 1 • f and p! — 1 • R, where a circumflex denotes 
the corresponding unit vector. The resonance condition between the two points is 
that the line separation At; is equal to the relative line-of-sight velocity along the 
ray, 

Aw = (v(r) - v(R)] ■ 1, (57) 

which, for radial outflow, becomes 

Aw = /zw(r) — p'v(R). (58) 

The resonance point radius R = R(r , p, Aw) may be fouml as the solution of equation 
(58), eliminating p' through the law of sines, R 2 (l - p' 2 ) = r 2 (l - p 2 ). 

We shall consider only the particular case of the velocity law (44) with p — 

In this case we find that R is a root of the cubic, 

R 3 - R.R 2 + (l - p 2 ) r 2 (R. - R) = 0. (59) 
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In solving this equation, we seek real rootB such that R > R.. For fixed r and Av, 
there are no such roots when — 1 < fx < n m { n , where 


Mmin ~ 


Aw - V c 


v{r) 


For all other jz, there is exactly one such root, but for n* < fx < Umax i where 


l^max — 


Av 

•w 


(60) 


(61) 


this root is associated with a resonance point that is obscured by the stellar core. 
When an unobscured resonance point does exist, then the intensity in the blue wing 
of the red line is given by equation (56); otherwise it is given by equation (13). 

The distance l that radiation travels between scatterings in the blue and red 
line can, in general, be on the order of a stellar radius or more. Since, in the local 
stability analysis conducted here, we are not concerned with effects that require 
perturbations to He coherent over such large distances, we ignore any effect that 
perturbations in the blue-line intensity might have on the red line. We therefore 
retain from §11 the blue wing boundary condition /i r (r, — » oo) = 0 for the 

perturbed red-line intensity. We do, however, modify the mean intensity to reflect 
these multiple scattering effects. 

For the particular case that the blue line is optically thick (assumption 3), we 
find from equations (16) and (23) that 


Iob(R,n\x— -oo) = S ob (R) =- 


(62) 


where 0b and 0 c b are evaluated by equations (24) and (25) for the blue line. Note 
that, although at fixed R, Sob is not a function of p,', in general we have R = 
i?(r,/z, Av), and so S 0 t> = -S'obC J2 (m))- Using equations (56) and (62), we may then 
generalize the function D(fi, r) (cf. equation (13)) to describe the overall angle 
dependence of the intensity in the blue wing of the red line, 


£> r (/i,r) = 


f°, 

0cb(R)/Pb(R ), 

i, 

,0cb(R)/0b(R), 


-1 < IX < IX min 
Mmin < fX < fX* 
(X* < IX < IXmax 
Mmai < fX < 1. 


(63) 


The second equality represents the effect of backscattering by the blue line, whereas 
the last equality represents the effect of shadowing and subsequent forward scat- 
tering of the stellar core radiation. With these two blue-line scattering effects thus 
incorporated into the definition of D, the remainder of the perturbation analysis 
then proceeds as in §11. We thereby find that the red-line perturbed force g t r in also 
given by equation (43) if we just replace D by D r and 0 e by 0 er = ( D r (fx)p(fx )) (cf. 
equation (25)). These two modifications represent generalizations of the “direct” 
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and “drag” terms of the single-line perturbed force to include the additional effects 
of scattering by the overlapping blue line. 

As we are primarily interested here in the effects of line overlap on the insta- 
bility to optically thin perturbations, we restrict our discussion to results for the 
perturbation growth rate 5t(<7i/vi) in the short wavelength limit k » x- For a 
wind with the velocity law given by equation (44) with p = 5 , figs. 4a and 4b 
show the radial variation of the red-line growth rate 9?(ffi r /vi) for line pairs with 
velocity-unit wavelength separations Av ranging from 0.1 Voo to 2 voo. For the case 
of maximum separation Av = 2vooi the lines no longer overlap, and so we recover 
in this case the earlier result (cf. figure 1) for an isolated single line. Furthermore, 
since the blue line of a pair is never influenced by the red line, this isolated line 
result also gives the perturbed blue-line growth rate 5R(^it/ui). 

For velocity separations in the range < Av < 2 ^ (figure 5b), note that the 
red-line growth rate actually increases with increasing line overlap. This is because, 
in this separation range, the dominant effect of line overlap is to enhance the “direct” 
term in equation (43) through backscattering by the blue line of radiation from the 
opposite side of the stellar wind. For Av < Uoo (figure 5a), on the other hand, 
this direct term is diminished beyond the radius at which shadowing by the blue 
line reduces the intensity from the stellar core; the red-line growth rate at large 
radii thus now decreases with increasing overlap. For strongly overlapped lines 
(An £ 0.2Voo)> the combined effects of reduction in the direct term by shadowing 
and enhancement in the drag term by backscattering causes the net red-line growth 
rate to become slightly negative very near (r ^ 1.1J2*) the stellar surface. Since 
the growth rates for unshadowed lines all become very small near the surface, this 
result implies that the combined effect of an ensemble of overlapping lines would 
tend to make the very base of the stellar wind stable, (See however the discussion 
at the end of §IVa.) However, only slightly away from the surface (r £ 1.1 fib), the 
growth rates for even the strongly shadowed lines (e.g., Av = 0.1v©o) again become 
positive, although they are still reduced relative to the growth rates for unshadowed 
lines. Since the growth rates for these unshadowed lines are all very large away from 
the wind base, the combined effects of an ensemble of shadowed and unshadowed 
lines is to make the bulk of the wind very unstable. 

We may thus summarize the results of our analysis of line overlap by stating 
that, while the contribution of strongly shadowed lines to the flow instability (as 
well as to the flow driving) is indeed reduced compared to unshadowed lines, their 
overall effect in the outer parts of the wind is not, as claimed by Lucy (1984), one 
of damping. Since the unshadowed lines (e.g., the bluest line of a multiplet) still 
make a strong contribution to the growth rate, we can conclude again that such 
radiatively driven flows are very unstable. 


V. SUMMARY AND CONCLUSIONS 

In this paper, we have extended our earlier study (paper I) of the stability of 
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line-driven stellar winds to take proper account of the finite angular size of the 
stellar photosphere and to include the dynamical effects of scattered radiation. Our 
work differs from the previous study of scattering by Lucy (1984) in that it is 
based on the physically more reasonable assumptions that the scattering occurs 
with complete redistribution (CRD) (vs. complete coherence) and that the mean 
wind flow is spherically symmetric (vb. plane parallel). Our results are summarized 
as follows: 

1) The mean scattered radiation exerts a drag force on velocity perturbations, 
which can eliminate the contribution of scattering lines to the instability of the 
flow (Lucy 1984), but only at the base of the wind, and only if one ignores 
photospheric limb darkening. With limb darkening, the scattering-line growth 
rate at the wind base is positive, and is a few percent of the rate for pure 
absorption lines (derived in paper I). 

2) As the flow moves outward from the base, the scattering-line growth rate rapidly 
increases to more than 50% of the rate for pure absorption lines within a stellar 
radius of the stellar surface, and eventually reaches 80% of that rate far from 
the star. 

3) Perturbations in the scattered radiation field change the details of the “bridging 
law”, which describes the response to intermediate wavelength perturbations, 
but this change is such as to affect mostly the propagation, and not the growth, 
of perturbations. 

4) The contribution of the red component of an overlapping line pair to the wind 
instability may be reduced compared to that of isolated lines, but the overall 
effect is not one of damping in the bulk of the wind. 

The basic conclusion we draw from these results is that the bulk of material in 
a radiatively driven stellar wind is still very unstable, even after all the effects of 
scattered radiation (including photon drag and line overlap) are taken into account. 


It U a pleasure to acknowledge helpful discussions with D. Friend and K. Mac- 
Gregor. This work was supported by the NASA HEAO-2 Guest Investigator Pro- 
gram and by NASA grant NSG-7406. 


APPENDIX: B(ijl) AND ASSOCIATED AVERAGES 

In order to evaluate the various angular averages occuring in equations (37) 
through (40), we must be able to evaluate the quantities A[n) and B[n), defined 
in equations (35) and (36). We first note that the double integral defining B(fi) is 
identical to that in the expression for Sg/Sv in paper I (equation |I-24]). In paper 
I a good approximation for 6g/6v was presented in equation (1-28), implying the 
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corresponding approximation 


BM « 


r M + 


Here the blue-edge absorption frequency x M is a Blowly-varying function of t m , given 
as the solution to the equation 

= 1. (A2) 

In the case of a Doppler profile for r M » 1, d>(x M ) « 2x ^/r^, leading to the approx- 
imation 

* r „|i*w( 2x.)r M3) 

Although is a function of p, we ignore this weak dependence in equation {A3} 
and take x M « x. = constant in evaluating the angular averages. 

We do not have a similar approximate expression for A(p), so we must find other 
methods to evaluate it. Other methods are also necessary if it is desired to obtain 
more accurate values of B(fi) than those given by equation (Al). A convenient way 
to evaluate both A(/x) and i?(/z), as well as their angular averages, is to express 
them in terms of Fourier integrals. 

Let us define = l/r M = /i 2 7i 4- (1 - p 2 )-] ro> in terms of which K ^ 
where K = k/x o is the wavenumber normalized by the mean line opacity xo- By 
changing the variable of integration to t = (x' - x)/7/i in the inner integrals in 
equations (35) and (36), and interchanging the orders of integration, we obtain 

A(m)= [°° dt e~ iK>lt K (*)7m)» M<) 

Jo 

B(/x) = r 7 „), (A5) 

Jo 

where 

/C(t,7) = J dx<f>{x)<f>(x + 7 *)exp {~7 -1 |${x) - $(* + '7*))} • i AQ ) 

£{ l n) = [ dx 4>(x)<f>(x + 7<) exp ( — 7 _1 ^(x)} . (A7) 

J — oo 

The various angular averages of A(p) and B{fx) can now easily be expressed as 
Fourier-type integrals. For example, 

(A(/i)> = § J-^dfiAifi) = \ [A(/z) + A(-/x)]. (A8) 

Since 7 M is even in /x, we have from equation (A4), 

(A) = [ d/x f dt cos K fxt JC(t,'f M ) 

Jo Jo 

r oo rl fin i 

— l dt cos Kt / — JC(— ,7u), (A9) 

Jo Jo u u 
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using the transformation t —* t/p and interchanging the order of integration. Simi- 
larly, 

(pA) — f Q dt sin Kt dp ( /V 10 ) 

From the values of (A) and (pA) the quantity r?, defined in equation (40), can be 
found. 


The required averages of B arc somewhat more complicated, since each involves 
both even and odd functions of p. These averages are most simply obtained using 
the approximations (Al) or (A3), but they can also be expressed in terms of Fourier- 
type integrals in a manner analogous to the derivation of equations (A9) and (A 10 ). 
We omit the details but give the results: 




ikp 2 
Q o 


B{p) 




= iK[ dte ,Kt h f dp—D[p)£(~ } ‘ 7 ^) 
Jo 1 J». ip p 

-il< [“dt cos Kt~ [ l dp -£(-,% M ), 
Jo (3 Jo ~fn (X 

= iK [ dte~ tKt ^ [ dp-D(p)£{-,'in) 
Jo *Jn. 'in M 

fOQ B Lt? t 

-Kj dt sin/ff ^ [ dp— £(-,>), 

Jo p Jo /i 


(/ 111 ) 


(A12) 


The representation of these angular averages by means of Fourier integrals has 
several advantages: l) The existence of standard numerical methods for evaluating 
Fourier integrals allows closer control of the cancellation of rapidly varying complex 
exponential or trigonometric terms, so that high accuracy can be obtained. 2) Since 
the wavenumber K appears only in the final Fourier transform, the numerical work 
of evaluating the two inner integrals is much reduced. 3) Complex arithmetic, if 
required at all, comes in only during the last step of evaluating the Fourier integrals. 
4) The asymptotic behavior at small and large wavenumbers can be found from 
well-known general results on Fourier integrals. 

Let us now find the asymptotic behavior of A and B for large and small values 
of K , using the representations (A4) and (A5). For K — ► 0 (Sobolev limit) we find 

roo 

A -* dt K (t, in) = 1 — P{d)i (A 13) 

Jo 

where p{p) is defined in equation (18). In obtaining equation (A13) we have used 
equation (A 6 ) and the fact that d$(z) = —<j>(x)dx. Similarly for K — > 0 we find, 

B — > p{p) — e" r **. (A14) 

In this limit A is purely real, so the imaginary part of (pA) vanishes. Thus »? — ► 0 > 
implying that the perturbed scattered intensity field vanishes in the Sobolev limit, 
as was expected on the basis of simple arguments. 
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For K -* oo (optically thin perturbation limit), equations (A4) and (A5) yield 
the results A —* 0 and B — ♦ 0, using the Riemann-Lcbcsguo lemma. This implies 
rj — > 0, so that the perturbed scattered intensity field also vanishes in this limit. 
However, in the angular averages of equation (42), B appears in the combination 
ilxKpB, which requires a more accurate estimation. We therefore integrate equation 
(A5) by parts once to obtain 

inKB = £( 0,^) - (A15) 

Applying the Riemann-Lebesgue lemma for K —* oo and using K ^ — Kt we find 

ipK^B -* r /(t £(0,7 M ). (A 16 ) 

We may also use the approximate equation (A3) to express this limit as 

inK^B^—. (A17) 
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FIGURE CAPTIONS 

Figure 1: Radial variation of the line growth rate (in units of the growth rate uj. 
defined in text) for short wavelength, optically thin perturbations (fc x)> as given 
by the perturbed line force to velocity ratio g\fv\. The upper (lower) curves are 
for a pure absorption (scattering) line, and the index p measures the steepness of 
the assumed velocity law. 

Figure 2: Wavenumber variation of the peturbed scattering ratio tj (defined by 
equation (40) in the text) for a line with radial optical depth t x = 100 at various 
radii r in a wind with a velocity law given by equation (44) with p = 

Figure 3: Wavenumber variation of real (solid curves) and imaginary parts of per- 
turbed force to velocity ratio for a scattering line with radial optical depth T\ — 100 
at various radii r/R t in the wind. Curves labeled So = Si =0 denote results for 
pure extinction (no scattering) case, while curves labeled Sj = 0 denote results if 
perturbations in an existing scattered field are ignored. 

Figure 4: Stability-dispersion curves showing the wavenumber variation of the real 
(solid) and imaginary (dashed) parts of the frequency for outward (+) and 
inward (— ) propagating radiative acoustic waves with various values of the sound 
speed y/ia. 

Figure 5: Same as figure 1, except for the red component of an overlapping pair of 
scattering lines for various values of the velocity separation Av. 
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